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Abstract 

Neutrino oscillations in matter can exhibit a specific resonance enhancement - 
parametric resonance, which is different from the MSW resonance. Recently it has 
been shown that the oscillations of atmospheric and solar neutrinos inside the earth 
can undergo parametric enhancement when neutrino trajectories cross the core of the 
earth. In this paper we continue the study of the parametric resonance of neutrino 
oscillations in the earth. The following issues are discussed: stability of the resonance 
with respect to the variations of the zenith angle of the neutrino source; higher-order 
resonances; prospects of the experimental observation of the parametric resonance of 
neutrino oscillations. We also comment on a recent controversy regarding the physical 
nature of the resonance enhancement of the oscillations of the core crossing neutrinos 
in the earth. 
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1 Introduction 



It is well known that neutrino oscillations in matter can differ significantly from oscillations 
in vacuum, the best studied example being the Mikheyev-Smirnov-Wolfenstein (MSW) effect 
H], |2|. The MSW effect is an enhancement of the neutrino flavor transition probability due 
to an amplification of the oscillation amplitude in matter. However, possible enhancement 
effects of matter on neutrino oscillations are not limited to the MSW effect. It has been 
pointed out in |3], |J] that the probability of neutrino flavor transition can also be strongly 
enhanced if the oscillation phase undergoes certain modification in matter. This can happen 
if the change of the matter density along the neutrino path is correlated in a certain way 
with the change of the oscillation phase. This amplification of the neutrino oscillation 
probability in matter due to specific phase relationships has an interesting property that 
it can accumulate if the matter density profile along the neutrino path repeats itself, i.e. 
is periodic. The phenomenon is analogous to the resonance in dynamical systems whose 
parameters periodically vary with time - parametric resonance. It was therefore called 
parametric resonance of neutrino oscillations |3], |J. In particular, in |4j] the parametric 
resonance was considered for neutrino oscillations in a matter with a periodic step function 
("castle wall") density profile, which allows an exact analytic solution. 

The parametric resonance can lead to large probabilities of neutrino flavor transition in 
matter even if the mixing angles both in vacuum and in matter are small. This happens 
because each half-wave oscillation of the transition probability is placed on the top of the 
previous one, i.e. the transition probability "drifts up". If mixing angle in matter is very 
small (matter density is far from the MSW resonance density), the parametric resonance 
enhancement of neutrino oscillations can manifest itself only if the neutrinos pass through 
a large number of periods of density modulation, i.e. travel a sufficiently long distance. 
However, if matter density is not very far from the MSW resonance one, an interesting 
interplay between the MSW effect and parametric effects can occur. In particular, a strong 
parametric enhancement of neutrino oscillations can take place even if the neutrinos pass 
only through 1-2 periods of density modulation ||. 

Although the parametric resonance of neutrino oscillations is certainly an interesting 
physical phenomenon, it requires that very special conditions be satisfied. Unfortunately, 
these conditions cannot be created in the laboratory because this would require either too 
long a baseline or neutrino propagation in a matter of too high a density (see Sec. 5 below). 
Until recently it was also unclear whether a natural object exists where these conditions 
can be satisfied for any known source of neutrinos. This situation has changed with a 
very important observation by Liu and Smirnov || (see also |7j]), who have shown that 
the parametric resonance conditions can be approximately satisfied for the oscillations of 
atmospheric into sterile neutrinos v s inside the earth. 

It is known that the earth consists of two main structures - the mantle and the core. 
Within the mantle and within the core the matter density changes rather slowly (the density 
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variation scale is large compared to the typical oscillation lengths of atmospheric and solar 
neutrinos), but at their border it jumps sharply by about a factor of two. Therefore to 
a good approximation one can consider the mantle and the core as structures of constant 
densities equal to the corresponding average densities (two- layer model). Neutrinos com- 
ing to the detector from the lower hemisphere at zenith angles in the range defined by 
cos = (—1) -T- (—0.837) traverse the earth's mantle, core and then again mantle. Therefore 
such neutrinos experience a periodic "castle wall" potential, and their oscillations can be 
parametrically enhanced. Even though the neutrinos pass only through "1.5 periods" of 
density modulations (this would be exactly one period and a half if the distances neutrinos 
travel in the mantle and in the core were equal), the parametric effects on neutrino oscil- 
lations in the earth can be quite strong. Subsequently, it has been pointed out by Petcov 
H that the parametric resonance conditions can also be satisfied (and to even a better 
accuracy) for the z/ 2 <-> v e oscillations in the earth in the case of the v e - v^( T ) mixing []. 
This, in particular, may have important implications for the solar neutrino problem. The 
parametric resonance in the oscillations of solar and atmospheric neutrinos in the earth was 



further explored in a number of papers ||, ID[ [11]. In the present paper we continue the 



study of this phenomenon. The following issues are addressed: the stability of the resonance 
with respect to the variation of the zenith angle of neutrino trajectory (Sec. 3), parametric 
resonances of higher order (Sec. 4), and the prospects of experimental observation of the 
parametric resonance in neutrino oscillations (Sec. 5). In the last section we comment on 
a recent controversy regarding the physical nature of the resonance enhancement of the 
oscillations of the core crossing neutrinos in the earth. 



2 Parametric resonance in neutrino oscillations in the 
earth 

In this section we shall summarize the main formulas describing neutrino oscillations in the 
earth which will be relevant for the subsequent discussion. We will be essentially following 
the formalism and notation of ref. || . 

Consider the oscillations in a 2-flavor neutrino system in a matter with periodic step 
function density profile. We will be assuming that one period of density modulation consists 
of two parts of the lengths R m and R c , with the corresponding effective matter densities N m 
and N c ("castle wall" density profile). For the v e <-> v^t) oscillations the effective matter 
density coincides with the electron number density, whereas for the v e ,n,T *-* v s oscillations 
in an isotopically symmetric matter it is a factor of two smaller. The parametric resonance 
in such a system occurs when the oscillations phases 20 m and 2<p c acquired over the intervals 

1 In H a different name for this phenomenon was suggested - neutrino oscillation length resonance. We 
prefer to follow the original terminology of H |j. 
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R m and R c are odd integer multiples of tt j|, ^ 0. Let us denote 



S = -^-, V i = ^ N i> u> i = y /(co S 2e 6-V i ) 2 + (sm26 6y (i = m, c) . (1) 

Here E, Am 2 and #o are the neutrino energy, mass squared difference and vacuum mixing 
angle, respectively. The difference of the neutrino eigenenergies in a matter of density iV; is 
2uJi, so that the oscillations phases acquired over the intervals R m and R c are 

20 m = 2cj m R m , 2(f) c = 2uj c R c . (2) 

The parametric resonance conditions can be written as || [|, ||, [], 0, §, [| 

1 = | + A;7r, 02 = ^ + ^, k,k' = 0,1,2,... (3) 

The evolution of a system of oscillating neutrinos is conveniently described by the evolution 
matrix U. In the case of the 2-flavor system U is a 2 x 2 unitary matrix which can be 
parametrized as 

[/ = Z - ierW Z 2 + W 2 = 1 . (4) 

Here a are the Pauli matrices in the flavor space. The evolution matrix U describes the 
evolution of an arbitrary initial state and therefore contains all the information relevant for 
neutrino oscillations. In particular, the probabilities of the neutrino flavor oscillations P 
and of u 2 <-» v e oscillations P 2e are given by 

P = W 2 + W 2 , P 2e = sin 2 fl + W 1 (W l cos 26 + W 3 sin 26 ) . (5) 



We shall consider neutrino oscillations in the earth in the two-layer approximation. 
Neutrinos coming to the detector from the lower hemisphere of the earth at zenith angles 
6 in the range cos 6 = (—1) (—0.837) (nadir angle n = 180° — < 33.17°) traverse the 
earth's mantle, core and then again mantle, i.e. three layers of constant density with the 
third layer being identical to the first one. In this case the parameters Z, W take a very 
simple form ||: 

Z = 2 cos0 m Y — cos 0c , (6) 

W = (2 sin m sin 29 m Y + sin C sin 26 c , 0, — (2 sin m cos 20 m Y + sin C cos 26* c )) . (7) 

Here 9 m and 9 C are the neutrino mixing angles in the matter of the density N m and N c (i.e. 
in the mantle and in the core of the earth), 

Y = cos m cos 0c — sin m sin C cos(26 l m — 29 c ) , (8) 

2 In refs. |], || these conditions were derived for the particular case k — k' , which, however, includes 
the most important principal resonance with k = k' = 0. 
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and the vector W was written in components. At the parametric resonance, i.e. when the 
conditions (|^) are satisfied, the neutrino flavor transition probability takes the value Jj, [7] , 

P = sm 2 (26 c -46 m ), (9) 

whereas the probability of the v-i <-> v e transitions is |§ 

P 2e = sin 2 (2# c - A6 m + 6 ) . (10) 

These probabilities can be close to unity (the arguments of the sines close to 7r/2) even 
if the amplitudes of neutrino oscillations in the mantle, sin 2 29 m , and in the core, sin 2 29 c , 
are rather small. This can happen if the neutrino energy lies in the range E c < E < E m , 
where E m and E c are the values of the energy that correspond to the MSW resonance in the 
mantle and in the core of the earth. The probability P 2e is relevant for the description of the 
oscillations of solar neutrinos in the earth fT2|, [13| . In the case of small mixing angle MSW 
solution of the solar neutrino problem, sin 2 2#o < 1CT 2 [14], and P 2e practically coincides 
with P unless both probabilities are very small. 



The trajectories of neutrinos traversing the earth are determined by their nadir angle 
n . The distances R m and R c that neutrinos travel in the mantle (each layer) and in the 
core are given by 



R m = R [cos 9„ - \Jr 2 /R 2 - sin 2 Q n j , R C = 2R \jr 2 /R 2 - sin 2 Q n . (11) 

Here R = 6371 km is the earth's radius and r = 3486 km is the radius of the core. The 
matter density in the mantle of the earth ranges from 2.7 g/cm 3 at the surface to 5.5 g/cm 3 
at the bottom, and that in the core ranges from 9.9 to 12.5 g/cm 3 (see, e.g., JT3J). The 
electron number fraction Y e is close to 1/2 both in the mantle and in the core. Taking the 
average matter densities in the mantle and core to be 4.5 and 11.5 g/cm 2 respectively, one 
finds for the u e v^ T oscillations involving only active neutrinos the following values of 
V m and V c : V m = 8.58 x 1CT 14 eV, V c = 2.19 x 1CT 13 eV. For transitions involving sterile 
neutrinos u e <-» v s and v^ T *-> v s , these parameters are a factor of two smaller. 



3 Stability of the parametric resonance conditions 

If the parametric resonance conditions (|3[) are satisfied, strong parametric enhancement of 
the oscillations of core crossing neutrinos in the earth can occur || [7| || |9], [10], [TTJ 0. We shall 
now discuss these conditions. The phases 4> m and C depend on the neutrino parameters 



3 For the parametric resonance to be a maximum of the transition probability, an additional condition 
has to be satisfied (two conditions in the case of v-i v e oscillations ||). For the oscillations between 
neutrinos of different flavor, this condition can be written as cos^c — AB m ) < 0. For small vacuum mixing 
angles 80 it essentially reduces to 5 < V c , with a small region around S — V m excluded. 
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Am 2 , 8q and E and also on the distances R m and R c that the neutrinos travel in the mantle 
and in the core. The path lengths R m and R c vary with the nadir angle; however, as can be 
seen from (|TT1) , their changes are correlated and they cannot take arbitrary values. Therefore 
if for some values of the neutrino parameters a value of the nadir angle n exists for which, for 
example, the first condition in eq. (|3|) is satisfied, it is not obvious if at the same value of 0„ 
the second condition will be satisfied as well. In other words, it is not clear if the parametric 
resonance conditions can be fulfilled for neutrino oscillations in the earth for at least one 
set of the neutrino parameters Am 2 , #0 and E. However, as was shown in [§, ||, not only 
the parametric resonance conditions are satisfied (or approximately satisfied) for a rather 
wide range of the nadir angles covering the earth's core, they are fulfilled for the ranges of 
neutrino parameters which are of interest for the neutrino oscillations solutions of the solar 
and atmospheric neutrino problems. In particular, the conditions for the principal resonance 
(k = k' = 0) are satisfied to a good accuracy for sin 2 28 < 0.1, S ~ (1.1 1.9) x 10~ 13 eV 2 , 
which includes the ranges relevant for the small mixing angle MSW solution of the solar 
neutrino problem and for the subdominant «-> v e and v e <-> v T oscillations of atmospheric 
neutrinos. 

The fact that the parametric resonance conditions can be satisfied so well for neutrino 
oscillations in the earth is rather surprising. It is a consequence of a number of remarkable 



numerical coincidences. It has been known for some time |1^, [], |17j that the potentials V m 
and V c corresponding to the matter densities in the mantle and core, the inverse radius of the 
earth R~ l , and typical values of S = Am 2 /AE of interest for solar and atmospheric neutrinos, 
are all of the same order of magnitude - (3 x 10~ 14 - 3 x 10~ 13 ) eV. It is this surprising 
coincidence that makes appreciable earth effects on the oscillations of solar and atmospheric 
neutrinos possible. However, for the parametric resonance to take place, a coincidence by 
an order of magnitude is not sufficient: the conditions ([3]) have to be satisfied at least within 
a 50% accuracy 0. This is exactly what takes place. In addition, in a wide range of the 
nadir angles B„, with changing 0„ the value of the parameter 5 at which the resonance 
conditions @ are satisfied slightly changes, but the fulfillment of these conditions is not 
destroyed. We shall show now that in this row of mysterious coincidences, at least the last 
one - the stability of the parametric resonance conditions with respect to variations of the 
nadir angle - has a simple explanation. 

We begin with simple qualitative arguments. Consider the principal parametric reso- 
nance for which k — k' — 0. In figs. 1 and 2 the flavor transition probability P and the 
parameters m — n/2 and (fi c — n/2 are plotted as functions of 5 for sin 2 29q = 0.01 and 
n = 27.3°. The highest peak of P corresponds to 5 ~ 1.55 x 10~ 13 eV at which, as can be 
seen in fig. 2, both m — tt/2 and <p c — tt/2 nearly cross the abscissa, i.e. the conditions (]3|) 
are approximately satisfied. This is the parametric resonance peak; the smaller peaks on its 
left and right are due to the MSW resonances in the mantle and in the core, respectively. 
The minima of the curves <p m — 71/ 2 and <p c — tt/2 correspond to the minima of uo m and 
uj c . The values of R m and R c determine the steepness of the curves. With increasing nadir 
angle, R m increases and R c decreases, which means that the branches of the <fr m — tt/2 curve 
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become steeper and those of the 4> c — tt/2 curve become less steep. Increasing the steepness 
of 4> m — tt/2 would push the intersection point of the <p m — tt/2 and 4> c — tt/2 curves upwards; 
however, decreasing the steepness of <ft c — tt/2 would push this intersection point in the 
opposite direction. As a result, in a wide range of the nadir angles n , with changing Q n 
the intersection point shifts along the abscissa but remains close to it, i.e. the parametric 
resonance conditions are not destroyed. It is easy to see that this only happens when the 
intersection takes place between the minima of <p m — tt/2 and <p c — tt/2, i.e. for 5 in the 
range 

V m / cos 26 < 5 < V m / cos 26 . (12) 

If the intersection takes place outside this interval, the changes of the steepness of the 
4>m — 71 "/2 and (f> c — tt/2 curves due to changing O n will push their intersection point in the 
same direction, and this will quickly destroy the parametric resonance conditions. Therefore 
outside the range (|12"D the parametric resonance can only occur in rather narrow intervals 
of O n . Interestingly, the interval flE^) is not only favored due to the wide range of allowed 
values of n , but also corresponds to the region where the strongest parametric enhancement 
of neutrino oscillations takes place f| |j. 

We shall now quantify the above arguments. Each of the two parametric resonance 
conditions in (H) can be solved with respect to the value 5 of the parameter 5 at the 



resonance: 



5 = cos29 V m ±J[(2k + l)-K] 2 /4R 2 m - S m 2 29 V* (13) 



5 = cos 20 o K± V[( 2A;/ + l )A 2 /^ R l - sin 2 2fl \/ c 2 (14) 

The parametric resonance can occur if the values of 5$ given by these two expressions coincide 
or nearly coincide. Notice that the existence of the radicals in ([13|) and fll4]) in general puts 
upper limits on the vacuum mixing angle 8 . In particular, for the principal resonance, 
assuming that the resonance condition can be satisfied for all neutrino trajectories that 
cross the core (including the vertical ones), from (0) one finds 



sin 2 20 o < Z7 ^^^O.O4. (15) 



If one excludes small nadir angles, the constraint becomes less stringent. For example, for 
sin 2 n > 0.12 one obtains sin 2 29 < 0.07. The analogous upper bound following from flT3] ) 
are less restrictive except for the trajectories close to the edge of the core of the earth. 

The (±) signs in eqs. (|13"D and ( |14"D correspond to different regimes of the parametric 
resonance. We shall now choose the plus sign in ( |T3"D and the minus sign in (|I~4"D, which for 
small #o corresponds to the interval (12). The other cases can be considered similarly 0. 



4 In the cases that correspond to the two plus signs, the parametric resonance does not lead to a maximum 
of the transition probability; for two minus signs there exist parametric maxima, but they are typically 
smaller than in the region of 5 that we are considering. 
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The consistency of eqs. (|13|) and flip) requires 



cos 2^0 (V c — V rn 



[{2k + 1)tt] 



-sin 2 2fl + 



[(2A;' + 1)tt 



2 



For small vacuum mixing angles 8q this simplifies to [] 



-sin 2 29 I/ 2 . (16) 



/t, ,n 7T /2fc + 1 2k' + l\ . . 

(V c -V m ) = -[ —— + ——- . (17) 

Eq. ( |16"D was obtained from the requirement that both conditions in (|3]) are satisfied at the 
same value of 5. If it is fulfilled, with changing G n the resonance value of 5 changes, but 



the parametric resonance conditions are still satisfied. Thus, it is the condition (16) that 
ensures the stability of the parametric resonance conditions with respect to the variations 
of the nadir angle. It remains now to check if or when this condition is indeed satisfied. 
We shall do that in the limit of the small 9q when the simplified form ([H]) applies. Similar 
analysis can be performed with the use of eq. flTB]). 

The l.h.s. of eq. ( |TTD is 9 n independent, so must be the r.h.s. As we have pointed out 
above, with increasing n the path length R m increases whereas R c decreases, which is just 
a consequence of the spherical geometry of the earth. One can therefore expect at least a 
partial cancellation of the changes of the (2k + 1)/R m and (2k' + 1)/R C terms in their sum 
on the r.h.s. of eq. (|17D . We shall show now that for the principal parametric resonance 
with k = k' = this cancellation is indeed quite good in a rather wide range of G n . In 
fig. 3 the nadir angle dependence of -R" 1 , R~ l and (-R" 1 + R^ 1 ) is plotted. One can see 
that while i?" 1 and R~ x change noticeably with Q n , their sum remains almost constant in 
the interval Q n = -f- 0.4. For larger values of O n this sum increases, but the increase does 
not exceed 25% for O n < 0.54. Thus, if the parametric resonance conditions are satisfied 
at some value of O n in the interval < O n ^ 0.54, they will be approximately satisfied in 
the whole interval. This explains the stability of the parametric resonance conditions. The 
fact that they are satisfied to a good accuracy for at least one value in the above-mentioned 
interval is a consequence of the numerical coincidences that we discussed earlier. 



4 Higher-order parametric resonances 

Eq. (|16|) can be used for a systematic search for the parametric resonances of neutrino 
oscillations in the earth. Once a solution is found (i.e. the values of R m and R c are 
determined), one can find the corresponding value of 5q from (|1^) or ([HI). The properties 
of the principal resonance have been discussed in some detail in || [7|, || ||; here we shall 
concentrate on the higher-order resonances. 

5 For the particular case k = k' = this condition was obtained in B. 
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The resonances of interest are located in or close to the region (|T2|) . The resonances 
that are far from that region either do not correspond to maxima or lead to small values of 
the transition probability. In addition, they can exist only in rather narrow intervals of the 
nadir angle. We have been able to identify two higher-order resonances of the v e <-> z/ M ( T ) 
oscillations. Let us discuss their properties and possible implications. 

(i) The resonance with k — 0, k' — 1. The resonance enhancement of the transition 
probability occurs in the interval of the nadir angles G n ~ 20° 4 25° and for relatively 
large vacuum mixing angles, sin 2 29 ~ 0.68 4 0.86. It corresponds to 5o ~ (7 4 9) x 10~ 14 
eV, depending on the values of ra and sin 2 2#o- For typical values of these parameters the 
resonance is illustrated by figs. 4 and 5. One can see that there is practically no detuning, 
i.e. the conditions @ are satisfied very well. The parametric peak at So — 9 x 10~ 14 eV is 
high and very wide, the energy width at half height being AE/Eq ~ 1 (the energies E\ and 
E 2 at the low-energy and high-energy borders of the half- height region satisfy E 2 /Ei > 3). 
The parametric resonance conditions can also be fulfilled for sin 2 2# < 0.68, but in this 
region the resonance corresponds to a low-lying saddle point rather than to a maximum of 
the transition probability. The existence of the resonance with k — 0, k' — 1 was pointed 
out in [PH| , but its properties were not discussed there. 

(ii) The resonance with k — 0, k' — 2. This resonance occurs for small values of the 
nadir angle, 0„ £ 12°, and large (close to maximal) vacuum mixing, sin 2 29q > 0.9. The 
resonance values of S lie in the range 5q ~ (0.65 4-1.1) x 10~ 13 eV. With increasing sin 2 29 
the resonance value of S decreases. This resonance is shown for Q n = and sin 2 29 = 0.95 in 
figs. 6 and 7. The parametric peak at <5 ~ 1 x 10~ 13 eV is slightly less wide than in the case 
(i): E2/E1 ~ 2.1. One can see from fig. 6 that the curves <p m — tt/2 and <f) c — 5ir/2 cross the 
abscissa at different but close points which are located symmetrically around 5 = 1 x 10~ 13 
eV. The detuning at 5 = 1 x 10~ 13 eV does not exceed ±0.15, i.e. is rather small. 

What are possible implications of these higher-order resonances? Unlike the principal 
resonance which can occur for both small and large values of #o; the resonances (i) and (ii) 
are only possible if the vacuum mixing angle is rather large, sin 2 29q > 0.7. For atmospheric 
neutrinos, for values of Am 2 of practical interest, such large mixing in the v e — v„ system 



is excluded by the Super-Kamiokande []18| and CHOOZ |L9|] data. Large mixing angle is 
in principle of interest for the vacuum oscillations and large mixing angle MSW solutions 
of the solar neutrino problem. However, the former requires very small values of S, orders 
of magnitude below those at which the parametric resonance of the neutrino oscillations 
in the earth can occur. The large mixing angle MSW solution, on the contrary, requires 
too large values of 5. For Am 2 ~ 7 x 10~ 6 eV 2 which is the smallest value allowed by 
the total event rates of the four solar neutrino experiments ||14|| , and E — 10 MeV which 
is a typical boron neutrino energy, one gets 5 ~ 1.75 x 10~ 13 eV. This is about a factor 
of 1.6 - 2.5 larger than the values of 5q that are necessary for the higher-order parametric 
resonances discussed above. These resonances could, however, occur for hep neutrinos which 
constitute the highest-energy component of the solar neutrino flux (endpoint energy 18.8 
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MeV). Smaller values of Am 2 (and, consequently, of 8) are allowed if the boron neutrino flux 
is considered as a free parameter |TJ, |2IJ. Relatively small values of Am 2 are also allowed 
in 3-flavor analyses of the solar neutrino data, even when the constraints on the mixing 



matrix element U e 3 following from the CHOOZ experiments are taken into account |2T^ [20 
Notice, however, that although the values Am 2 ~ 7 x 10 -6 eV 2 or even less are allowed 
at large vacuum mixing angles by the total event rates, they are disfavored by the recent 
Super-Kamiokande data on the zenith angle distribution of the solar neutrino events. 

We conclude that the higher-order parametric resonances are less interesting than the 
principal one. The domain of the neutrino parameters in which they can occur for the 
oscillations of solar neutrino in the earth is currently disfavored by the solar neutrino data, 
though not completely ruled out. These resonances are of no interest for the atmospheric 
neutrino oscillations. It would be interesting to look for other possible neutrino sources for 
which these resonances could play a role. 



5 Can the parametric resonance of neutrino oscilla- 
tions be observed? 

Besides being an interesting physical phenomenon, the parametric resonance of neutrino 
oscillations can provide us with an important additional information about neutrino prop- 
erties. Therefore experimental observation of this effect would be of considerable interest. 
We shall now discuss the prospects for experimental observation of the parametric resonance 
of neutrino oscillations in the earth, having in mind mainly the principal resonance. There 
are two main sources of neutrinos for which the parametric resonance can be important - 
atmospheric neutrinos and solar neutrinos. Both sources have their advantages and disad- 
vantages from the point of view of the possibility of observation of the parametric resonance. 
We shall now briefly discuss them. 

We start with atmospheric neutrinos. The parametric resonance can occur in the <->• 
v s |], [7j and also in the subdominant v e <-> channels of oscillations [10|. It can affect 



the distributions of /i-like events and also (in the case of the v e <-> v^ir) oscillations) lead to 
interesting peculiarities in the zenith angle distributions of the multi-GeV e-like events. 

The observation of the parametric effects is hampered by the loose correlation between 
the directions of the momenta of atmospheric neutrinos and of the charged leptons which 
they produce and which are actually detected. Because of this the trajectories of neutrinos 
coming to the detector are not known very precisely. In addition, the data are presented for 
certain samples of events (sub-GeV, multi-GeV, upward through-going, upward stopping) 
which includes collecting data over rather wide energy intervals. The contributions of the 
parametric peaks may therefore be integrated over together with other possible enhancement 
peaks - due to the MSW resonances in the mantle and in the core, making the distinction 
between these effects difficult. Also, strong resonance enhancement effects (both parametric 
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and MSW) can only occur either for neutrinos or for antineutrinos, depending on the sign of 
Am 2 []. The present atmospheric neutrino experiments do not distinguish between neutrinos 
and antineutrinos, therefore possible matter effects are "diluted" in the sum of the v- and 
^-induced events. At certain values of the ratio of the muon and electron neutrino fluxes 
r(E, Q n ) depending on the value of mixing angle #23 the parametric effects on e-like events 
are suppressed jlOfl . 

Atmospheric neutrinos have some advantages for observation of the parametric reso- 
nance of neutrino oscillations in the earth. Neutrinos come to the detectors from all di- 
rections, which means that practically the whole solid angle covering the earth's core will 
contribute to the effect. There are no additional suppression factors due to a specific com- 
position of the incoming neutrino flux which may quench the earth's matter effect on the 
oscillations of solar neutrinos (see below). Parametric effects may provide a unique probe 
of the neutrino mixing angle #13 with sensitivity going far beyond that of the long-baseline 



accelerator and reactor experiments EDI |22|. Possible ways of improving the prospects of 
the experimental observation of the parametric effects in the atmospheric neutrino oscil- 
lations include using various energy cuts, finer zenith angle binning and detectors capable 
of detecting the recoil nucleon, which would enable one to reconstruct the direction of an 
incoming neutrino [ ID , 2"2" |. 



Solar neutrinos can experience a strong parametric enhancement of their oscillations in 
the earth if the small mixing angle MSW effect is the correct explanation of the the solar 
neutrino deficit f|, ||. The parametric enhancement can occur in a wide range of values of 
sin 2 2#o and for the nadir angles n almost completely covering the core of the earth. The 
trajectory of each detected neutrino is exactly known. For boron neutrinos the resonance 
occurs at the values of Am 2 which correspond to the central part of the allowed interval for 
the small mixing angle MSW solution of the solar neutrino problem. 

However, there are some disadvantages, too. Unfortunately, due to their geographical 
location, the existing solar neutrino detectors have a relatively low time during which solar 
neutrinos pass through the core of the earth to reach the detector every calendar year. 
The Super-Kamiokande detector has a largest fractional core coverage time equal to 7%. 
In [|TjJ it was suggested to build a new detector close to the equator in order to increase 



the sensitivity to the earth regeneration effect; this would also maximize the parametric 
resonance effects in oscillations of solar neutrinos in the earth. In the case of the MSW 
solutions of the solar neutrino problem the probability P$e of finding a solar v e after it 
traverses the earth depends sensitively on the average u e survival probability in the sun Ps 
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1 2P 

p S e = p s + — ~ sin2 °o) ■ ( 18 ) 

cos 2u 

The probability P 2e can experience a strong parametric enhancement, but in the case of 



6 In general, the parametric enhancement can take place even for "wrong sign" Am 2 , but in the case of 
neutrino oscillations in the earth these effects are small. 
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small mixing angle MSW solution of the solar neutrino problem the probability Ps for 
the Super-Kamiokande and SNO experiments turns out to be rather close to 1/2. This 
means that the effects of passage through the earth on solar neutrinos should be strongly 
suppressed. The current best fit of the solar neutrino data is not far from the line in the 
parameter space where Ps is exactly equal to 1/2 and Pse = Ps the earth matter 

effects are absent). Whether or not it will be possible to observe the parametric resonance 
of the oscillations of solar neutrinos in the earth depends on how close to this line the true 
values of sin 2 29 and Am 2 are. By now the Super-Kamiokande experiment has not observed, 
within its experimental accuracy, any enhancement of neutrino signal for earth core crossing 
neutrinos ||23|| . This can be because the parametric enhancement of the neutrino oscillations 



in the earth does not occur {e.g. if the true solution of the solar neutrino problem is vacuum 
oscillations or large mixing angle MSW effect), or because the values of Am 2 and sin 2 2# 
are too close to those at which Ps = 1/2. Hopefully, with accumulated statistics of the 
Super-Kamiokande and forthcoming data from the SNO experiment the situation will soon 
be clarified. 

It is interesting to note that the Super-Kamiokande data on the zenith angle dependence 
of the solar neutrino events seems to indicate some deficiency of the events due to the core- 
crossing neutrinos rather than an excess [23|, although it is not statistically significant. 



Should this deficiency be confirmed by future data with better statistics, it could have 
a natural explanation in terms of the parametric resonance of neutrino oscillations. As 
follows from (0), the parametric enhancement of P^ e for core crossing neutrinos can lead 
to a deficiency of the events if the neutrino parameters are in the small-sin 2 2#o part of 
the allowed region which corresponds to Ps > 1/2 (see, e.g., fig. 10 in ref. |24}| ). In this 
case one should also have an "opposite sign" overall day-night effect (fewer events during 
the night than during the day). In any case, given the current experimental constraints 
on the neutrino parameters, if the small mixing angle MSW effect is the true solution of 
the solar neutrino problem, the only hope to observe earth matter (day-night) effects on 
solar neutrinos seems to be through the parametric resonance of oscillations of core crossing 
neutrinos. 

As we have seen, observing the parametric resonance in oscillations of solar and atmo- 
spheric neutrinos in the earth is not an easy task. Can one create the necessary matter 
density profile and observe the parametric resonance in neutrino oscillations in the labora- 
tory (i.e. short-baseline) experiments? Unfortunately, the answer to this question seems to 
be negative. The parametric resonance can occur when the mean oscillation length in mat- 
ter approximately coincides with the matter density modulation length 0, |], [5]: l m — L. 
In a matter of density Ni the oscillation length is given by l m = 7r/u>i where u>i was de- 
fined in ([I]). Let us require l m ^ 1 km. Assume first that Vi <; 5, i.e. u>i are dominated 
by matter density terms. Than for l m < 1 km one would need a matter of mass density 
Pi > 3.3 x 10 4 g/cm 3 , clearly not a feasible value. Conversely, for pi < 10 g/cm 3 , one finds 
l m > 3300 km, a distance comparable with the earth's radius. Consider now the opposite 
case, 5 ^> Vi. Then the oscillation length in matter essentially coincides with the vacuum 



11 



oscillations length which in principle can be rather short provided that the vacuum mixing 
angle 8q is small (otherwise this would contradict reactor and accelerator data). However, 
in this case there is another problem. Requiring l m < 1 km one finds 5 > 2.5 x 1CT 10 eV. 
For pi < 10 g/cm 3 one therefore has Vi/5 < 1CT 3 . This means that the mixing angles in 
matter are very close to the vacuum one, 6{ ~ 6q(1 + Vi/5), and so their difference is very 
small: AO = 6 2 — Q\ — (AV/S) 9$ < 10 _3 #o- When the difference of mixing angles in matter 
is small, the parametric effects can manifest themselves only if neutrinos travel over a large 
number of periods, n ~ tt/AA6 0, |], || p|. Therefore in this case the necessary baseline 
is ~ it 2 / (AOqAV) > 3 x 10 3 km, again too large. One can conclude that the sole presently 
known object where the parametric resonance of neutrino oscillations can take place is the 
earth, as was first pointed out in || [?J . 



6 Parametric resonance or oscillation length resonance? 

Recently, there have been a considerable controversy regarding the physical nature of the 
resonance enhancement of the oscillations of the core crossing neutrinos in the earth. In 
refs. H [7|, ||, |26|, [TOj this enhancement was interpreted as a parametric resonance effect. 



However, in refs. P, |27|, BR ITTJ it was interpreted as a new effect, called neutrino oscillation 



length resonance. It was argued that this effects has nothing to do with the parametric 
resonance of neutrino oscillations. This controversy has led to a considerable confusion and 
we therefore believe that it is worthwhile to discuss the issue in some detail. We shall argue 
here that the resonance enhancement of oscillations of core crossing neutrinos is a typical 
parametric resonance effect. 

The general possibility of the parametric resonance in neutrino oscillations was sug- 
gested independently in || and f| and further studied in || . In these papers it was shown 
that the parametric resonance is a specific resonance enhancement of neutrino transition 
probability which occurs due to a certain relationship between the density modulation scale 
and neutrino oscillation length, does not require large amplitude of neutrino oscillations, and 
has a property to accumulate when the matter density profile repeats itself, i.e. is periodic. 
The resonance enhancement of the oscillations of core crossing neutrinos in the earth fully 
matches this description. The conditions (11) in M (reproduced also in P7L E3, O, |29|), 



which are claimed to be the conditions for the oscillation length resonance, coincide with 
the parametric resonance conditions given in || [|, |5], || [7| (see eqs. (10) and (17) in j|, eqs. 
(2) and (4) in ||), and in particular eq. (7) in || and eq. (12) in @ where the paramet- 
ric resonance conditions were formulated specifically for the case of the oscillations of core 
crossing neutrinos in the earth). Moreover, the probability of the neutrino flavor transition 
at the oscillation length resonance given in |jTT| coincides with the oscillation probability at 
the parametric resonance given in & |7| (see eq. (Bh above). 



We shall now comment on the arguments presented in refs. ||, ^7], |28|, pj§ (and given 



below in italics) in support of the point of view that the resonance enhancement of the 
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oscillations of core crossing neutrinos in the earth is very different from the parametric 
resonance effects in the case of the "castle wall" matter density profile studied in 

• In /[J/ the results were presented for the specific case of small densities (V mfi <C 5) which 
does not correspond to the situation with oscillations of neutrinos in the earth. 

The case V m ,c "C 8 was discussed in ref. J| just in order to stress that the proposed effect 
- parametric resonance in neutrino oscillations - is very different from the MSW effect. The 
parametric resonance can take place for any relation between the matter-induced potentials 
and S (see, e.g., 0). 

• In ^ the case of the density profile with equal spatial dimensions of the layers of constant 
density (R m = R c ) was studied, whereas in the earth R m ^ R c . 

The nature of the effect does not depend on whether R m and R c are equal to each 
other or not. In fact, there is a value of the nadir angle (0„ = 27.3°) which determines 
the neutrino trajectory through the earth along which R m = R c . Fig. 1 was drawn for 
this particular value of the nadir angle; it shows a rather typical pattern of the resonance 
enhancement in the oscillations of core crossing neutrinos (compare with analogous figures 
in refs. 

• The density profile along the trajectory of core crossing neutrinos is not periodic even in 
the two-layer approximation; neutrino path length falls short even of 1.5 periods. 

We shall first comment on the second part of this argument. In general, it is incorrect: 
neutrino path length is less than 1.5 periods of density modulation for cosG n > 0.89 but 
exceeds this value for 0.837 < cos6 n < 0.89, which constitutes about 31% of the solid angle 
covering the earth's core. This, however, is unimportant. The term "1.5 periods" was used 
in refs. 0, |K|] just as a shorthand for describing the situation when neutrinos cross three 



layers of constant density with the third layer being identical to the first one. 

In fig. 8 we show the coordinate dependence of the probability of neutrino oscillations 
in the earth for B ra = 11.5° (solid curve), and for comparison also show the probability 
in the hypothetical case when neutrinos propagate through the full two periods of density 
modulation (dashed curve). Obviously, in both cases the physical nature of the enhancement 
of the transition probability is the same: the probability drifts up and each new half-wave 
piece of the oscillating curve is placed on the top of the previous one. This pattern is typical 
of the parametric resonance of neutrino oscillations |3|, |4|, ||. 

Now about the periodicity of the matter density profile of the earth in the two-layer 
approximation. In the strict mathematical sense, a function f(x) can be called periodic 
only if it is defined in the infinite interval x G (— oo, oo). It was shown, however, in |3], |J 
that the complete flavor conversion is possible even if neutrinos travel over a finite distance 



7 Similar figures can also be found in some earlier papers |30| | the authors of which have found the 
parametric resonance of neutrino oscillations in the earth in numerical calculations but have not recognized 
its parametric nature. 
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in a matter of periodic density. In || it was demonstrated that this distance can be quite 
short (1-2-3 periods) if the parametric resonance energy is not very far from the MSW 
resonance energy. This is exactly what happens in the case of the oscillations of the core 
crossing neutrinos in the earth. The term "periodic" was used in refs. p, |4], |5], [| in this 
limited sense - as a finite piece of the periodic density profile. Since the path length of 
core crossing neutrinos in the earth exceeds one period of density modulation R m + R c , the 
density profile of the earth in the two-layer approximation satisfies this definition. 

Actually, the only relevant question here is if it is really important that the third layer 
traversed by the neutrinos going through the earth is identical to the first one. This brings 
us to the next objection. 

• The specific enhancement of the neutrino transition probability under discussion can take 
place even if neutrinos cross a number of layers of constant but different densities. 

In general, this is correct. Assume, for example, that the neutrinos cross three layers of 
densities N± < N 2 < N 3 , and that we are free to choose the spatial extensions of these layers, 
so that the condition similar to @ can be satisfied. Then, in principle, one can have a strong 
enhancement of the transition probability. However, the resonance conditions are just 
the extremality conditions; for the resonance to be a maximum of the transition probability, 
additional conditions have to be satisfied. In the case of neutrino flavor oscillations in 
a step-function density profile with = Aq, there is only one such condition; in the 
case Ni < N2 < A3 there are three additional conditions. Thus, in this latter case there 
are altogether six conditions (three conditions on oscillation phases and three additional 
conditions), whereas in the case A3 = Aq there are only three conditions. 

As was discussed above, the possibility of the resonance enhancement of the oscillations 
of the core crossing neutrinos in the earth is a result of a number of remarkable numerical 
coincidences. It is because of these coincidences that all three conditions for the parametric 
resonance of neutrino oscillations can be satisfied. It is much less probable that six different 
conditions could be satisfied as a result of some coincidences of parameters. Fortunately, in 
the case of the oscillations of the core crossing neutrinos in the earth, the third layer that 
the neutrinos cross is identical to the first one, i.e. the matter density profile is periodic 
and so there is a chance that the resonance enhancement of neutrino oscillations under the 
discussion will be observed. 

This work was supported by Fundagao para a Ciencia e a Tecnologia through the grant 
PRAXIS XXI/BCC/16414/98 and also in part by the TMR network grant ERBFMRX- 
CT960090 of the European Union. 
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Figure captions 



Fig. I. Transition probability P for v e <-> oscillations in the earth for core-crossing 
neutrinos vs 5. sin 2 29 = 0.01, 6 n = 27.3°. 

Fig. 2. Phases <j) rn — n/2 (dashed curve) and <p c — it/ 2 (solid curve) vs 5 for u e <-> u^ T 
oscillations in the earth, sin 2 29 = 0.01, n = 27.3°. 

Fig. 3. Parameters l/R m (long-dashed curve), 1/R C (short-dashed curve) and 1/R m + 1/R c 
(solid curve) as functions of the nadir angle 6 n . All distances are in the units of the earth's 
radius R; Q n is in radians. 

Fig. 4. Same as fig. 1 but for sin 2 29 = 0.77, 6„ = 23.3°. 

Fig. 5. Phases m — n/2 (solid curve) and C — 37r/2 (dashed curve) vs 5 for v e <-> v^ T 
oscillations in the earth, sin 2 29 = 0.77, n = 23.3°. 

Fig. 6. Same as fig. 1 but for sin 2 26 = 0.95, Q n = 0. 

Fig. 7. Phases <p m — ir/2 (solid curve) and 4> c — 57r/2 (dashed curve) vs 5 for u e <-> i/ M;T 
oscillations in the earth, sin 2 29 = 0.95, n = 0. 

Fig. 8. Solid curve: transition probability P for u e <-> z/ M T oscillations in the earth as a 
function of the distance t (measured in units of the earth's radius R) along the neutrino 
trajectory. 5 — 1.8 x 10~ 13 eV, sin 2 29 = 0.01, n = 11.5°. Dashed curve: the same 
for a hypothetical case of neutrino propagation over full two periods of density modulation 

(tmax = 2(-R m + Rc))- 
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